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NOMENCLATURE 


A  area,  m2 

Ac  area  of  finite  difference  grid  element  face  perpendicular  to 

cylinder  axis,  m2 

Apin  area  of  finite  difference  grid  element  face  parallel  to  cylinder 
axis  and  closest  to  center  of  cylinder,  m2 
Apout  area  of  finite  difference  grid  element  face  parallel  to  cylinder 
axis  and  farthest  from  center  of  cylinder,  m2 
a  thermal  diffusivity,  n^/s 

Cp  specific  heat,  J/g-°K 

c^  specific  heat  (liquid)  J/g-°K 

Cg  specific  heat  (solid)  J/g-°K 

dh  hydraulic  diameter,  2(R^  —  R0) 

E  internal  energy,  J/g 

Fo  Fourier  number,  ai/R2 

h  heat  transfer  coefficient,  W/m2-°K 

i  enthalpy,  J/g 

i*  enthalpy  in  the  solid  saturated  state,  J/g 

s 

k  thermal  conductivity,  W/m-°K 

thermal  conductivity  (liquid),  W/m-°K 
kg  thermal  conductivity  (solid),  W/m-°K 

ky  tube  wall  thermal  conductivity,  W/m-°K 

L  length  of  PCM  cylinder,  m 

Lj  latent  heat  of  fusion,  J/g 

M  number  of  radial  divisions  in  finite  difference  grid 

•  . 

m  flow  rate  of  liquid  metal,  g/s 


N 


number  of  axial  divisions  in  finite  difference  grid 


NPIPES 

Nux 

PCM 

Pe 

Pr 

q 

r 

R(t) 

Re 

Rir 


O 

Ste 

T 

T 

Td 

Tavg 


in 


init 


out 


number  of  PCM  cylinders 
Nusselt  number,  hx/k 
phase  change  material 
Peclet  number,  VL/a 
Prandtl  number,  u*Cp/k 
thermal  energy  flow,  J/s 
radial  coordinate,  m 

radial  position  of  the  phase  change  front,  m 

Reynolds  number,  vx/v 

radius  of  liquid  metal  boundary,  m 

density 

radius  of  PCM  cylinder,  m 

Stefan  number,  cg  (Tm  -  dimensionless 

temperature,  °K 

time,  seconds 

test  duration,  seconds 

average  of  maximum  inlet  and  outlet  liquid  metal  temperatures, 
°K 

maximum  inlet  temperature  of  liquid  metal,  ‘’K 
initial  temperature,  °K 

maximum  outlet  temperature  of  liquid  metal,  °K 
melting  temperature  of  PCM,  °K 
wall  thickness,  m 

Kirchoff  temperature,  defined  in  Sect.  3.1 
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ABSTRACT 


A  space-based  thermal  storage  application  must  accept  large 
quantities  of  heat  in  a  short  period  of  time  at  an  elevated  temper¬ 
ature.  An  approximate  model  of  a  lithium  hydride  phase  change  energy 
storage  system  was  used  to  estimate  reasonable  physical  dimensions  for 
this  application  which  included  the  use  of  a  liquid  metal  heat 
transfer  fluid.  A  finite  difference  computer  code  was  developed  and 
used  to  evaluate  three  methods  of  enhancing  heat  transfer  in  the  PCM 
energy  storage  system.  None  of  these  three  methods,  inserting  thin 
fins,  reticulated  nickel,  or  liquid  lithium,  improved  the  system  per¬ 
formance  by  a  significantly  large  amount.  The  use  of  a  95%  void  frac¬ 
tion  reticulated  nickel  insert  was  found  to  increase  the  storage 
capacity  of  the  system  slightly  with  a  small  decrease  in  the  system 
energy  density.  The  addition  of  10%  liquid  lithium  was  found  to  cause 
minor  increases  in  both  storage  density  and  storage  capacity  with  the 
added  benefit  of  reducing  the  hydrogen  pressure  of  the  lithium 


1 .  INTRODUCTION 


1 . 1  Project  Purpose 

A  space-based  power  system  application  requires  a  large  amount  of 
heat  rejection  from  a  power  cycle  in  a  short  period  of  time.  One  pro¬ 
posed  solution  is  to  use  a  phase  change  thermal  energy  storage  system 
to  absorb  the  rejected  heat  and  then  release  it  to  space  over  a  longer 
period  of  time.  The  basic  arrangement  of  the  proposed  storage  system 
consists  of  a  cylindrical  tube  filled  with  a  phase  change  material 
with  a  high  melting  temperature  surrounded  by  an  annular  region  con¬ 
taining  the  liquid  metal  heat  transport  fluid.  Lithium  hydride  is  the 
preferred  thermal  energy  storage  medium  because  of  its  high  latent 
heat  and  thermal  conductivity  and  because  its  melt  temperature  matches 
the  needs  of  the  application  well. 

Similar  systems  have  been  modeled  in  the  past.  However,  the 
short  time  constraints  for  this  application,  coupled  with  the  poor 
thermal  diffusivity  of  the  medium,  have  led  to  the  consideration  of 
three  heat  transfer  enhancements.  The  enhancements  considered 
include:  1)  the  addition  of  thin  fins,  2)  the  insertion  of  a  solid 

lattice  of  reticulated  nickel,  and  3)  the  addition  of  a  molten  metal 
to  a  like  metal  salt  (lithium  and  lithium  hydride).  In  this  third 
enhanced  system,  the  metal  is  always  a  liquid  and  the  salt  melts  and 
freezes,  forming  a  slush-like  material.  The  reticulated  nickel  lat¬ 
tice  is  similar  to  a  steel  wool  pad  inserted  into  the  PCM  tube  before 
the  PCM  material  is  poured  in. 
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This  research  is  part  of  a  larger  effort  with  other  researchers 
investigating  two  related  areas.  The  first  is  the  formation  of  voids 
within  the  phase  change  medium  caused  by  changes  in  density  (assumed 
constant  in  the  present  study)  and  the  influence  of  a  gravity-free 
environment  on  the  location  of  such  voids.  The  second  is  the  per¬ 
formance  of  a  packed  bed  of  encapsulated  phase  change  material. 

The  primary  focus  of  this  project  is  to  define  a  baseline  latent 
thermal  energy  storage  system  and  compare  its  performance  to  three 
enhanced  latent  thermal  energy  storage  systems.  An  existing  computer 
code  was  substantially  modified  to  reflect  the  requirements  of  this 
task.  Chapter  2  describes  an  approximate  analysis  that  was  used  to 
estimate  feasible  system  configurations.  In  Chap.  3,  the  computer 
model  of  the  system  is  described  and  the  results  for  the  enhanced 
systems  are  compared  in  Chap.  4.  The  conclusions  from  this  comparison 
and  recommendations  for  future  study  are  presented  in  Chap.  5. 

1 .2  Problem  Definition 

At  this  time,  the  energy  storage  requirements  for  the  application 
have  been  only  tentatively  defined  and  are  subject  to  change.  The 
power  source  may  range  in  size  from  10  to  5,000  MW  and  may  operate 
over  a  time  period  of  200  to  1,500  seconds  [1].  For  this  study,  a 
thermal  storage  system  was  sought  to  accept  the  rejected  heat  from  a 
constant  power  source  of  250  MWe  (generated  with  a  thermal  efficiency 
of  25%)  for  an  operational  period  between  900  and  1000  seconds.  The 
thermal  energy  that  must  be  stored  is  therefore  approximately 


7.5  x  1011  J.  The  storage  system  is  composed  of  a  system  of  tubes 
filled  with  PCM  that  receives  the  rejected  heat  from  hot  liquid  metal 
flowing  on  the  shell  side,  see  Fig.  1.1.  The  maximum  tube  length  per¬ 
mitted  is  about  5  m.  The  maximum  allowable  liquid  metal  temperature 
leaving  the  PCM  storage  system  is  1100°K  and  the  maximum  inlet  liquid 
metal  temperature  is  1200°K.  The  initial  temperature  of  the  PCM  sys¬ 
tem  is  700°K.  The  tube  walls  are  made  of  304SS  with  a  tube  wall 
thickness  of  about  1.27E-4  m  (0.005  in.)  [2].  The  properties  of  the 
liquid  metal,  tube  wall,  and  PCM  are  shown  in  Table  1.1. 

An  individual  PCM  cylinder  surrounded  by  liquid  metal  is  modeled 
for  this  analysis,  see  Fig.  1.2.  The  PCM  cylinders  are  radially  sym¬ 
metric  with  insulated  ends.  The  outer  boundary  of  the  liquid  metal 
annular  volume  is  assumed  to  be  adiabatic  by  symmetry.  Dwyer  pre¬ 
sented  a  full  discussion  on  the  use  of  such  an  equivalent  circular 
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Fig.  1.1.  System  Schematic 


Table  1.1*  Material 
Properties 


Property 

Value 

Liquid  metala 

k 

28.9  W/m-°K 

cp 

1.043  J/g-°K 

P 

7.401  x  10s  g/o3 

U 

0.161  g/m-s 

Lithium  hydride^ 

P* 

5.5  x  105  g/m3 

k  (solid) 

4.18  W/m°K 

k  (liquid) 

2.51  W/m°K 

Cp  (solid) 

8.33  J/g°K 

c  (liquid) 

7.36  J/g°K 

4 

2845  J/g 

T 

m 

962  °K 

Tube  wall‘d 

P 

7.803  x  106  g/m3 

k 

24.2  W/m-°K 

CP 

0.46  J/g-°K 

aNa  (567.)  K  (447.)  at 
704°C  (977°K)  [3,  p.  564]. 

b[4]. 


^Assumed  constant  for 

both  solid  and  liquid  states. 

^304SS  steel  (18%  C, 

8%  Ni)  (3,  p.  541]. 


boundary  for  triangular  arrays  of  tubes  with  wide  spacing  (Rf/RQ  > 
1.35)  [5,  p.  104-105].  His  methods  could  be  used  to  facilitate  the 
extension  of  the  results  of  this  study  to  a  full  bank  of  such  tubes. 
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1.3  Literature  Review 

The  use  of  a  hollow  cylinder  of  PCM  for  latent  energy  storage  was 
modeled  by  Solomon  et  al.  and  is  described  in  Ref.  6.  He  used  a 
finite  difference  formulation  with  the  Kirchoff  temperature  to  cal¬ 
culate  the  internal  energy,  temperature,  and  the  position  of  the  phase 
change  front.  The  Kirchoff  temperature  approach  is  also  used  in  this 
analysis  and  is  described  more  fully  in  Sect.  3.1.  Solomon  also  used 
a  quasi-steady  approximation  method  to  estimate  reasonable  values  for 
input  to  the  more  exact  finite  difference  model.  A  very  similar 
quasi-steady  approximation  was  used  in  this  effort  and  is  reported  in 
Sect.  2.1.  Solomon’s  work  and  this  analysis  differ  in  four  major 
areas:  1)  geometry,  2)  heat  transfer  material,  3)  the  charging 

requirements,  and  4)  the  use  of  heat  transfer  enhancements.  Solomon's 
geometry  consisted  of  a  hollow  cylinder  of  PCM  with  a  gaseous  heat 
transfer  medium  flowing  through  the  center  of  the  cylinder.  The 
charge  and  discharge  cycles  were  attuned  to  a  solar  cycle.  Because 


the  charging  times  were  so  long,  the  poor  conductivity  of  the  PCM  was 
of  less  concern  and  enhancements  were  not  considered. 

The  enthalpy  method  of  analysis  was  rigorously  verified  for  a 
phase  change  application  by  Shamsundar  and  Sparrow  [7].  Using  an 
integral  formulation,  they  show  that  the  enthalpy  approach  is  identi¬ 
cal  to  the  more  classical  approach  and  that  knowledge  of  the  location 
of  the  phase  change  front  is  not  necessary  for  the  solution. 

Shamsundar  and  Sparrow  also  defined  dimensionless  parameters  useful 
for  phase  change  problems  and  formulated  the  conservation  equation  in 
terms  of  these  variables.  These  variables  are  described  in  Sect.  3.1 
and  are  used  in  the  presentation  of  results  in  Sect.  4.2. 

A  useful  text  in  the  area  of  finite  difference  problem  solutions 
was  written  by  Dusinberre  [8].  This  text  includes  an  especially  clear 
discussion  of  time-step  size  requirements  for  stability. 

Knowles  and  Webb  [9]  have  examined  the  use  of  a  mixed  enthalpy 
model  for  metal/PCM  composities  both  analytically  and  experimentally. 

Their  thermal  model  is  based  on  volume-average  effective  thermal 
parameters  which  are  used  in  Sect.  3.4  of  this  report. 

Several  heat  transfer  correlations  for  liquid  metal  flowing  in  an 
annulus  were  available  in  the  literature.  Some  of  these  were  based  on 
theoretical  considerations  and  others  on  experimental  data. 

Dwyer's  theoretical  derivation  for  laminar  flow  was  based  on  the 
following  assumptions:  heat  transfer  from  the  inner  wall  only,  uni¬ 
form  heat  flux,  negligible  axial  heat  conduction,  and  fully  developed 
velocity  and  temperature  profiles  [10].  This  analysis  showed  the 
Nusselt  number  and  bulk  fluid  temperature  to  be  a  function  of  the 
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ratio,  Rf/R0‘  The  results  were  presented  in  tabular  form,  some  of 
which  are  listed  in  Table  1.2  (excerpted  from  Dwyer,  [10],  p.  337). 
The  simplistic  relation  between  geometry  and  Nusselt  number  in  this 
correlation  made  it  an  attractive  choice  for  the  approximate 
analysis.  These  results  were  therefore  used  in  the  approximate 
analysis  in  the  form  of  Eq.  2.10  in  Sect.  2.2  for  Rj/Rq  ratios  from 
1.2  to  3.0. 


Table  1.2.  Nusselt  Numbers  and  Wall 
and  Bulk  Temperature  Relationships 
for  Flow  Through  Annuli 


V'o 

Nu 

t>  -  V 

Tl  -  T2 

1  (parallel  plates) 

5.385 

0.743 

1.25 

5.60 

0.758 

1.50 

5.78 

0.771 

2.00 

6.17 

0.790 

4.00 

7.78 

0.832 

6.00 

9.22 

0.856 

8.00 

10.57 

0.870 

10.00 

11.90 

0.879 

00  (round  pipes) 

OO 

1.000 

Excerpted  from  Dwyer  [10, 
p.  337]  for  laminar  flow  with  heat 
transferred  from  the  inner  wall  only. 


Ti  -  inner  wall  temperature, 
T^  -  liquid  metal  bulk  temperature, 
and  T2  “  outer  wall  temperature. 


Dwyer  ([10],  p.  340)  also  provided  a  correlation  (Eq.  1.1)  appli¬ 
cable  to  either  laminar  or  turbulent  flow  in  an  annulus  under  condi¬ 
tions  of  fully  developed  flow  and  uniform  heat  flux. 


Nu  -  4.82  +  0.697  y  +  0.0222  (tf/Pe) 


(0.758  y0*053) 


(Eq.  1.1) 
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&S 


where 


y  -  Rf/R, 


iji  «  dimensionless  average  value  of  ratio,  eH/em 


eu  “  K  /pc  «  eddy  diffusivity  o.  heat  transfer 

H  e  p 


e  “  u./p  “  eddy  diffusivity  of  momentum  transfer 
m  e 


Kg  *  eddy  thermal  conductivity 


Ug  -  eddy  dynamic  viscosity 


Another  heat  transfer  correlation  was  based  on  empirical  data 


(Duchatelle  [11,12,13])  for  a  heat  exchanger  with  NaK  flowing  in  an 


annular  space.  The  reported  experimental  results  cover  a  range  of 


Peclet  numbers  from  500  to  3000.  The  authors  extrapolate  the  relation 


down  to  a  Peclet  number  of  40  in  comparing  their  results  to  those  of 


other  authors.  This  correlation  was  chosen  for  use  in  the  finite 


difference  code  because  this  range  is  appropriate  for  the  modeled  flow 


with  a  Peclet  number  of  about  200.  In  Duchatelle's  correlation 


(Eq.  1.2),  the  Nusselt  and  Peclet  numbers  are  based  on  the  hydraulic 


diameter  [2(R£-RQ)]  and  the  liquid  metal  physical  properties  are  taken 


at  the  bulk  average  temperature.  Duchatelle  also  uses  a  temperature 


dependent  conductivity,  unlike  the  constant  value  used  in  this 


study.  The  value  F  in  Eq.  1.2  is  a  factor  to  increase  the  heat 


transfer  coefficient  in  the  thermal  and  hydraulic  establishment  length 


is  shown  in  Eq.  1.3.  For  an  Rf/RQ  ratio  of  2  and  a  value  of  1, 


Eqs.  1.1  and  1.2  are  nearly  equivalent. 


Nu  -  (6.15  +  0.02  Pe0*8)  F  for  40  <  Pe  <  3000 


(Eq.  1.2) 
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(Eq.  1.3) 


F  -  1 


Another  equation  (Eq.  1.4)  was  theoretically  derived  for  liquid 
metal  flow  across  a  flat  plate  but  showed  good  comparison  to  experi¬ 
mental  data  for  a  flow  of  liquid  metal  along  the  surface  of  an  un¬ 


broken  section  of  tubing  [14). 


Nu  -  0.46  (Pe)0,65  for  103  <  Pe  <  105 


(Eq.  1.4) 


Yet  another  correlation  was  presented  by  Subbotin  [15,16]  for  a 
forced-convection  sodium  film  coefficient  for  a  bank  of  tubes, 

Eq.  1.5.  In  this  formula,  the  variable  "S"  is  the  distance  from  the 
center  of  one  tube  to  the  center  of  an  adjacent  tube.  The  ratio  of  S 
to  2Rq  is  therefore  analogous  to  the  ratio  Rf/R0  in  the  model  studied 
in  this  report.  The  correlation  in  Eq .  1.5  is  applicable  for  1.1  < 
S/2Rq  <  1.4  and  for  200  <  Pe  <  1200.  There  are  also  some  restrictions 
on  the  length  to  diameter  ratio.  For  the  geometry  of  the  configura¬ 
tions  modeled  in  this  report,  the  Subbotin  Nusselt  number  estimate  is 
about  25%  less  than  that  based  on  Duchatelle's  correlation.  However, 
the  Rf/R0  ratio  in  this  configuration  equals  three  and  is  beyond 
Subbotin's  limits.  If  configurations  with  smaller  radius  ratios  are 
modeled,  this  correlation  should  be  considered. 


Nu  ■  8 


\jr+ °-27  (s;- 


(Eq.  1.5) 


2.  PRELIMINARY  CONFIGURATION  SELECTION 


An  approximate  analytical  solution  for  this  problem  was  used  to 
determine  potential  system  configurations.  The  parameters  that  must 
be  coordinated  in  such  configurations  include:  PCM  cylinder  radius, 
liquid  metal  annulus  radius,  liquid  metal  velocity,  tube  length,  and 
number  of  cylinders  required. 

The  PCM  and  liquid  metal  regions  are  estimated  separately  in 
Sects.  2.1  and  2.2  and  the  two  approximate  solutions  are  then  coupled 
in  Sect.  2.3  through  their  common  geometrical  boundaries  and  the  con¬ 
vective  heat  transfer  coefficient.  The  approximate  solution  resulted 
in  a  small  self-contained  computer  code,  independent  of  the  finite 
difference  code.  It  should  be  noted  that  the  performance  of  the 
resulting  configurations  was  evaluated  using  the  finite  difference 
code.  Modifications  to  the  configurations  were  then  made  based  on 
this  evaluation. 

2. 1  Approximate  Model  for  the  Phase  Change  Material  Region 

The  PCM  approximate  model  is  patterned  very  closely  after  the  one 
introduced  by  Solomon  [6J.  Solomon  used  a  quasi-steady  formulation  to 
allcM  a  simplified  analysis.  This  approach  is  discussed  more  fully  in 
Arpaci  and  Larsen  [17,  p.  181].  A  quasi-steady  formulation  involves 
neglecting  unsteady  terms  in  at  least  one,  but  not  all,  of  the  govern¬ 
ing  equations  of  an  unsteady  problem.  The  application  of  the  quasi¬ 
steady  simplification  is  most  appropriate  when  the  phase  growth  is 
very  slow  and  the  sensible  heat  storage  is  much  smaller  than  the 
latent  heat  storage.  This  is  therefore  a  reasonable  approximation  for 


cases  where  the  Stefan  number  Is  less  than  or  equal  to  one  (see  Arpaci 
and  Larsen  [17,  p.  176]).  This  simplification  ignores  the  heat  addi¬ 
tion  required  to  steadily  heat  the  liquefied  layer  as  it  grows  in 
thickness. 

The  model  considered  consists  of  an  axisymmetric  PCM  cylinder  of 
outer  radius  RQ,  initially  solid  at  its  melting  temperature.  The 
surface  is  exposed  to  convective  heat  transfer  from  a  liquid  metal 
flowing  at  a  constant  temperature.  The  heat  transfer  coefficient  is 
assumed  to  be  a  known  quantity  during  this  portion  of  the  analysis  and 
is  discussed  more  fully  in  Sect.  2.3.  A  phase  change  front,  R(t), 
separates  the  liquid  PCM  (r  >  R(t))  from  the  solid  PCM  (r  <  R(t))  (see 
Fig.  2.1).  In  both  this  approximate  solution  and  in  the  finite  dif¬ 
ference  solution,  the  PCM  density  is  assumed  to  remain  constant,  i.e. 
the  density  is  unaffected  by  state  or  temperature.  Also,  there  is 
assumed  to  be  no  motion  (i.e.  no  convection)  within  the  liquid  PCM 
region.  Equations  2.1  through  2.8  follow  Solomon's  [6]  analysis 
except  for  changes  due  to  geometry. 
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Fig.  2.1.  Approximate  Model  of  the  PCM  Region. 
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The  temperature,  T(r,t),  In  the  liquid  region  satisfies  the  heat 
equation,  2.1. 


3T 

3 


T  a  3  /  3T\  „  v  . 

for  '  >  r(t) 


(Eq.  2.1) 


At  the  solid-liquid  interface,  the  temperature  equals  the  melting  tem¬ 
perature  and  the  rate  of  heat  melting  the  solid  PCM  (and  therefore 
moving  the  front)  equals  the  rate  of  heat  passing  through  the  front. 
Equation  2.3  is  discontinuous  at  R(t)  =  0  because  at  this  point  the 
front  ceases  to  exist. 


T(R(t) ,  x)  =  T 


m 


P  K 


3  R(t) 


3t 


-k 


3T(R(t),  t) 


3r 


for  R(t)  >  0 


(Eq.  2.2) 
(Eq.  2.3) 


At  the  container  wall  the  rate  of  heat  transfer  to  the  PCM  equals  the 
convective  heat  transfer  from  the  liquid  metal.  The  heat  capacity  of 
the  wall  is  assumed  to  be  negligible  here  and  is  discussed  further  in 
Sect.  3.3.2. 

3T(R  .  t) 


-k 


3r 


h(T(R  ,  t)  -  T  ) 
o  avg 


(Eq.  2.4) 


For  the  quasi-steady  approximation,  Eq .  2.1  is  replaced  by  the  steady 
state  equation  (Eq.  2.5).  It  should  be  noted  that  the  boundary  condi¬ 
tions  remain  time  dependent.  This  equation  must  be  constrained  to 
R( t)  >0  because  the  temperature  of  the  center  of  the  PCM  is  a  func- 


3T(0  t) 

tion  of  time  only  (i.e.,  ■  —  “  0)  and  the  quasi-steady  approxima¬ 


tion  fails  at  this  point. 


h  (‘  0-° 


for  r  >  R(t)  >  0 


(Eq.  2.5) 


The  solution  to  Eqs.  2.5  and  2.2  yields  Eq.  2.6. 


T(r ,  t)  -  T  +  a( x)  In  for  r  >  R(t)  >  0  (Eq.  2.6) 

m 

The  variable,  a(t),  is  found  by  substituting  this  solution  into 
Eq.  2.4  and  is  shown  in  Eq.  2.7. 


h(T  -  T  ) 
avg  m 


k  .  .  o 

R-  +  h  in 

o 


for  R(t)  >  0 


(Eq.  2.7) 


The  flux  equation  at  the  interface,  Eq .  2.3,  can  now  be  used  with 
Eq.  2.6  (recognizing  that  as  r  goes  from  0  to  t,  r  goes  from  Rq  to  R) 
to  solve  for  the  position  of  the  interface  as  given  implicitly  in 
Eq.  2.8. 


Up  to  this  point,  the  approximate  solution  in  the  PCM  region  is 
very  similar  to  that  shown  by  Solomon  [6],  From  here  onwards,  the 
solution  approach  will  be  different  to  reflect  the  differences  in  both 
the  problem  geometry  and  the  expected  process  application.  When  the 
PCM  cylinder  is  completely  melted,  R  is  equal  to  zero  (that  is,  the 
phase  change  front  has  moved  to  the  center  of  the  cylinder).  However, 
this  equation  is  discontinuous  at  R  =  0.  Therefore,  to  estimate  the 
radius,  R„,  for  which  the  PCM  cylinder  is  completely  melted  during  the 
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test  period,  R  is  set  equal  to  RQ/100  and  t  is  set  to  (test  dura¬ 
tion)  in  Eq.  2.8.  Solving  for  RQ  gives  Eq .  2.9. 


2  Ai, 

+  <T  ~  T  )  T. 

p  L  avg  m  d 


(Eq.  2.9) 


This  Rq  is  the  radius  that  will  yield  the  greatest  storage  dens¬ 
ity  for  the  given  set  of  simplified  assumptions.  For  tubes  with 
smaller  radii,  the  melting  is  also  complete  within  the  test  duration 
but  more  material  could  be  melted.  Therefore,  for  a  cylinder  of  fixed 
wall  thickness,  the  parasitic  mass  decreases  with  increasing  radius 
and  the  storage  density  can  increase  with  increasing  radius.  Above 
Rq,  some  of  the  PCM  present  will  not  melt  and  the  storage  density  will 
decrease. 

The  heat  transfer  coefficient  is  not  yet  known  and  must  be  esti¬ 
mated  before  Rq  can  be  calculated.  This  leads  to  the  analysis  of  the 
liquid  metal  region  and  the  coupling  of  the  two  regions. 


2 . 2  Approximate  Solution  for  the  Liquid  Metal  Region 

The  approximate  solution  for  the  liquid  metal  region  is  sought  to 
provide  guidance  in  the  selection  of  reasonable  physical  parameters 
for  the  finite  difference  solution.  As  discussed  in  the  last  section, 
an  estimate  for  the  convective  heat  transfer  coefficient  is  required 
to  estimate  the  PCM  cylinder  radius.  Also,  the  outer  radius  of  the 
liquid  metal  region  annulus  and  the  fluid  velocity  must  be  estimated 
to  avoid  exceeding  the  maximum  permissible  outlet  temperature. 

The  flow  of  liquid  metal  in  an  annular  region  was  examined  by 
Dwyer  [5,  p.  78-79  and  10,  p.  337-338]  for  a  set  of  simplifying 
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assumptions.  These  assumptions  include  1)  heat  transfer  from  the 
inner  wall  only,  2)  uniform  heat  flux,  3)  negligible  axial  heat  con¬ 
duction,  and  4)  fully  developed  velocity  and  temperature  profiles. 

This  problem  was  solved  exactly  for  the  laminar  flow  case  and  the  re¬ 
sults  for  average  Nusselt  number  and  bulk  temperature  are  presented  in 
tabular  form  (see  Sect.  1.3)  as  a  function  of  the  ratio,  R^/Ro.  These 
results  were  used  in  the  form  of  Eqs.  2.10  and  2.11  for  the  computer 
program  discussed  in  Sect.  2.3. 


Nu  =  0.78  Rc/R  +  4.6 
f  o 


h  = 


k  Nu 


(Eq.  2.10) 
(Eq.  2.11) 


The  bulk  liquid  metal  temperature  is  then  a  function  of  only  the 
wall  temperature,  the  thermal  flux,  and  the  radius  ratio,  Rf/RQ.  The 
wall  temperature  at  the  outlet  was  assumed  to  equal  the  PCM  melting 
temperature.  The  constant  heat  flux  was  taken  equal  to  the  required 
thermal  load  on  the  PCM  tube,  q.  This  gives  a  constant  liquid  metal 
temperature  for  a  constant  wall  temperature,  T  ,  as  shown  in 


Eq .  2.12. 


T.  =  T  +  q/(hA  ) 
Jtm  m  pout 


(Eq.  2.12) 


2. 3  Coupling  the  Liquid  Metal  and  Phase 
Change  Material  Approximate  Solutions 

The  PCM  approximate  model  defined  a  maximum  radius  for  which  the 
PCM  material  would  be  completely  melted  during  the  test  duration  as  a 
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function  of  the  heat  transfer  coefficient,  h.  The  liquid  metal  lami¬ 
nar  flow  solution  defines  the  temperature  and  heat  transfer  coeffi¬ 
cient  as  functions  of  the  tube  radius  and  the  ratio  of  the  adiabatic 
liquid  metal  radius  to  the  tube  radius.  These  two  solutions  are  ob¬ 
viously  tightly  coupled. 

Several  constraints  must  be  met  by  any  coupled  solution.  The 
liquid  metal  flow  and  temperature  change  must  meet,  without  exceeding 
greatly,  the  thermal  load  requirements  (Eq.  2.13).  The  mass  of  the 
PCM  and  the  liquid  metal  inventory  must  be  capable  of  storing  all  the 
rejected  heat  over  the  test  period  duration  (Eq.  2.14).  The  number  of 
PCM  tubes  required  should  be  as  low  as  possible.  The  liquid  metal 
velocity  was  constrained  to  values  less  than  5  m/s  because  higher  flow 
velocities  led  to  higher  liquid  metal  exit  temperatures. 


m  c(T  -  T(z-L))  >  W 
p  in  t 


(Eq.  2.13) 


[*(R  2  -  R  2)  L  (T  -  T  )  c  + 
f  o  out  init  p 

(tt  R  2  L)  (L  +  (T  -  T  )  c  )]  (NPIPES )  >  W  x,  (Eq.  2.14) 
o  t  m  inlt  p  t  d 

These  relationships  and  constraints  were  programmed  into  a  small 
computer  code  that  tested  a  wide  range  of  variables  to  define  feasible 
system  parameters,  listed  in  Table  2.1.  The  flow  chart  for  this  com¬ 
puter  code  is  shewn  in  Fig.  2.2. 
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Table  2.1.  Preliminary  Systen 
Parameters  Based  on 
Approximate  Analysis 


(n? 

L 

f  m) 

V 

(n/s ) 

NPIPF.S 

0.048 

0.124 

4.5 

0.05 

5000 

0.049 

0.136 

4.5 

0.04 

4600 

0.049 

0.136 

4.5 

0.09 

4600 

0.048 

0.124 

5.0 

0.10 

4  500 

0.049 

0.136 

5.0 

0.09 

4100 

0.044 

0.124 

5.0 

0.05 

5000 

0.053 

0.149 

3.5 

0.04 

5000 

0.053 

0.149 

3.5 

0.08 

5000 

0.053 

0.137 

4.0 

0.04 

4700 

0.053 

0.137 

4.0 

0.09 

4  700 

0.053 

0. 149 

4.0 

0.04 

4300 

0.053 

0.149 

4.0 

0.08 

4300 

0.053 

0.125 

4.5 

0.10 

4600 

0.053 

0.137 

4.5 

0.09 

4200 

0.053 

0.149 

4.5 

0.08 

3900 

0.053 

0.125 

5.0 

0.10 

4100 

0.053 

0.137 

5.0 

0.09 

3800 

0.053 

0. 149 

5.0 

0.08 

3500 

0.053 

0.149 

5.0 

0.12 

3500 

3.  FINITE  DIFFERENCE  SOLUTION 


3 . 1  Background 

The  approximate  solutions  included  several  simplifications  that 
do  not  reflect  the  physical  situation  such  as:  1)  no  axial  varia¬ 
tions,  2)  uniform  heat  flux,  3)  a  quasi-steady  state,  4)  a  constant 
inlet  liquid  metal  temperature,  and  5)  an  initial  state  of  solid  PCM 
at  its  melting  temperature.  An  analytical  approach  to  the  overall 
problem  would  be  too  complex  to  allow  an  exact  solution.  Therefore,  a 
numerical  approach  was  used  to  model  the  behavior  of  a  phase  change 
storage  system  charged  by  a  liquid-metal  flow.  Within  the  PCM  region, 
this  model  is  very  similar  to  that  used  by  Solomon,  et  al.  [6].  How¬ 
ever,  the  geometry,  element  definitions,  and  use  pattern  are  dif¬ 
ferent. 

The  PCM  cylinder  was  assumed  to  be  radially  symmetric  with  insu¬ 
lated  ends.  The  liquid  metal  flow  has  an  artificial  adiabatic  ex¬ 
ternal  boundary  as  discussed  in  Sect.  1.2.  The  liquid  metal  is 
assumed  to  exchange  heat  only  within  the  heat  exchanger  and  with  the 
thermal  load  source.  That  is,  outside  the  PCM  heat  exchanger,  the 
fluid  travels  in  an  insulated  pipe  so  that  there  is  no  upstream  or 
downstream  conduction.  An  initial  temperature  for  the  system  is 
specified. 

An  energy  balance  is  performed  for  each  element  in  the  finite 
difference  grid  described  in  Sect.  3.2.  The  resulting  equations  are 
solved  explicitly  within  the  PCM  region  and  implicitly  within  the 
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IV. 


liquid  metal  region.  Within  the  PCM  region,  the  equations  are  iden¬ 
tical  to  those  used  by  Solomon,  et  al.  [6]  and  are  based  on  energy 
conservation,  Fourier's  law,  and  the  equation  of  state,  shown  in 
Eqs.  3. 1-3.3.  The  equation  of  state  for  constant  specific  heats  for 
each  phase  is  shown  in  Eq .  3.4  (where  the  internal  energy,  E,  of  the 
solid  state  PCM  at  its  melting  temperature  has  been  set  to  zero). 


p  |i  +  7.q  -  0 


q  =  -k(T)  VT 


C(T) 


dE(T) 

dT 

+  —  E  <  0  (solid  phase) 
me  r 

s 


(Eq.  3.1) 
(Eq.  3.2) 
(Eq.  3.3) 


T  -J 


,  0  <  E  <  (melting  zone) 


(E  -  L  ) 

yj'm  + - c -  ,  E  >  (liquid  phase) 


(Eq.  3.4) 


A  "Kirchoff"  temperature,  U,  was  defined  by  Solomon  [6,  p.  18] 
and  is  useful  because  the  position  of  the  phase  change  front  becomes  a 
part  of  the  solution  rather  than  an  assumption  that  must  be  iterated 
upon  within  the  calculations.  Equations  3.5-3.10  outline  these 
energy  conservation  equations  in  terms  of  this  "Kirchoff"  temperature 
(defined  in  Eq.  3.6). 


K(T) 


k  ,  T  <  T 
s  m 


kL  *  T  >  Tm 


(Eq.  3.5) 
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(Eq.  3.6) 


VU  -  K(T)VT 


(Eq.  3.7) 


q  =  -VU 


(Eq.  3.8) 


P  ||  -  V.(VU) 


(Eq.  3.9) 


,  E  <  0 


0  ,  0  <  E  <  L 


k£  (E-Lf) 


»  E  >  L 


(Eq.  3.10) 


This  Kirchoff  temperature  approach  used  by  Solomon  is  very  simi¬ 


lar  to  the  enthalpy  model  discussed  by  Shamsundar  and  Sparrow  [7], 


They  rigorously  verified  the  enthalpy  model  and  defined  a  dimension¬ 


less  enthalpy,  0,  as  shown  in  Eq.  3.11.  The  dimensionless  time  param¬ 


eter  is  the  Fourier  number,  Fo,  Eq .  3.12.  Shamsunder  and  Sparrow  used 


the  dimensionless  product,  FoSte  as  a  dimensionless  time  parameter 


because  it  correlates  the  results  for  various  values  of  Ste,  defined 


in  Eq.  3.13.  The  enhanced  PCM  mixtures  discussed  in  Chap.  4  represent 


a  range  of  Stefan  numbers,  so  this  product  has  been  used  in  the  pre¬ 


sentation  of  results. 


1  r 

pV  h  P 


(i  -  i  *) 


(Eq.  3.11) 


21 


Fo  -  crr/R  2  (Eq.  3.12) 

o 

Ste  -  c  (T  -  T.  .„)/L,  (Eq.  3.13) 

s  m  inlt  t 


For  a  material  with  an  assumed  constant  density  and  a  constant 
latent  heat,  and  because  the  internal  energy  variable,  E,  was  defined 
equal  to  zero  at  the  solid  state  melting  temperature,  Eq.  3.11  can  be 
rewritten  as  Eq.  3.14.  In  this  equation,  9  describes  the  energy  con¬ 
tent  of  the  PCM  relative  to  the  latent  storage  capacity  and  relative 
to  the  energy  content  at  the  solid  state  melting  temperature.  A 
dimensionless  measure  of  the  energy  storage  in  the  PCM  can  be  con¬ 
structed  by  taking  the  difference  between  the  value  of  9  at  a  given 
point  in  time  and  the  initial  value  of  9,  see  Eq.  3.15.  The  dimen¬ 
sionless  temperature,  <(>,  is  defined  in  Eq .  3.16.  The  results  will  be 
presented  in  terms  of  these  dimensionless  variables  for  generality. 


l  1  M’N 

9  ■  VL^  A, EdV  ■  VE7  ^  Vii 


9*  =  9  -  9 


initial 


T  -  T 


init 


T  —  T 
m  init 


(Eq.  3.14) 
(Eq.  3.15) 
(Eq.  3.16) 


Within  the  liquid  metal  region,  no  phase  change  occurs.  Hence, 
the  Kirchoff  temperature  offers  no  advantage  and  the  internal  energy 
is  calculated  implicitly.  The  liquid  metal  calculations  were  ori¬ 
ginally  done  in  an  explicit  manner  similar  to  those  in  the  PCM 
region.  However,  the  liquid  metal  has  a  much  higher  conductivity  than 
the  PCM  and  required  a  very  small  time  step  to  avoid  stability  prob¬ 
lems  (see  Sect.  3.3.1).  This,  in  turn,  caused  the  computer  program  to 


consume  excessive  CPU  time  for  each  test  case.  This  portion  of  the 
code  was  therefore  rewritten  using  an  implicit  solution  method  for 
tridiagonal  matrices  outlined  in  Forsythe  and  Wasow  [18,  p.  104]. 

3.2  Finite  Difference  Grid  Definition  and  Solution 

The  problem  was  defined  in  two  spatial  coordinates,  r  and  z,  and 
one  temporal  coordinate,  t.  Although  each  node  represents  a  small, 
thin,  cylindrical  volume,  the  finite  difference  formulation  is  rec¬ 
tangular  in  nature.  The  grid  used  is  shown  in  Fig.  3.1.  The  areas 
and  volumes  of  these  grid  elements  are  shown  in  Table  3.1.  In  this 
table,  "i"  is  used  to  denote  the  node  index  in  the  r-direction  and  "j" 
is  used  for  the  z-direction.  A  single  r-direction  grid  node  is  used 
for  the  liquid  metal  region  because  this  is  sufficient  to  determine 
the  bulk  fluid  temperature  as  a  function  of  z  and  r.  This  bulk  fluid 
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Fig.  3.1.  Finite  Difference  Element  Grid. 


Table  3.1.  Finite  Difference  Areas  and  Volumes 
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f  o  2 
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2TTR  A£ 
o  2 

0 

TT(R  2~R  2)-^ 
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aArea  of  element  face  perpendicular  to  the  Z-axis. 

^Area  of  element  face  parallel  to  the  Z-axis  and  closest  to  the  center  of 
the  cylinder. 

^Area  of  element  face  parallel  to  the  Z-axis  and  farthest  from  the  center 
of  the  cylinder. 

^"M+l"  represents  the  node  in  the  liquid  metal  region. 

Area  taken  as  zero  because  this  is  an  adiabatic  surface. 


temperature,  in  turn,  is  sufficient  to  determine  the  heat  transfer  to 
the  PCM  cylinder  nodes  and  to  test  for  the  maximum  allowable  liquid 
metal  exit  temperature. 

The  internal  energy,  temperature,  and  Kirchoff  temperature  of 
each  node  are  taken  to  represent  both  the  point  values  and  the  mean 
values  for  each  node  at  each  point  in  time. 

A  PCM  node  is  shown  in  Fig.  3.2.  The  general  form  of  the  energy 
balance  is  shown  in  Eq.  3.17,  where  the  volumes  and  areas  are 
defined  in  Table  3.1.  The  superscript  of  "n+1"  indicates  the  value  at 
the  next  point  in  time.  Values  without  superscripts  indicate  the 
value  at  the  current  point  in  time.  The  updated  value  is  explicitly 
calculated  based  on  the  values  for  the  current  point  in  time.  The  use 
of  this  general  form  equation  is  defined  in  Table  3.2  where  qi ,  q2 , 
q 3 ,  and  q4  are  defined  in  terms  of  qa,  q^,  qc>  q^,  and  qg  (defined  in 
Eqs.  3.18-3.22)  for  all  nodes  within  the  PCM  region.  Equation  3.23 
demonstrates  the  use  of  the  general  form  (Eq.  3.17)  and  Table  3.2  for 
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Fig.  3.2.  PCM  Node  Energy  Balance. 
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+  a7  (ui,j+i-uij)  (E«‘  3*23) 
The  liquid  metal  region  is  solved  implicitly  using  Eq.  3.24  (for 
1  <  j  <  N)  and  3.25  (for  j=N).  The  areas  and  volumes  for  these  equa¬ 
tions  were  shown  in  Table  3.1.  After  this  matrix  of  equations  is 
solved,  the  liquid  metal  inlet  temperature  is  calculated  using 
Eq.  3.26. 
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(Eq.  3.26) 


T  (M+1,1)  =  T  (M+l.N)  +  T=- 

mcp 

The  heat  transfer  coefficient  between  the  liquid  metal  flow  and 
the  PCM  is  calculated  conservatively  (the  value  is  actually  greater  in 
the  thermal  and  hydraulic  development  length)  using  Eq .  3.27  [11]. 

This  and  other  correlations  available  for  liquid  metal  annular  flow 
cases  were  discussed  in  Sect.  1.3. 

Nu  =  (6.15  +  0.02  Pe°* 8)  (Eq.  3.27) 


3. 3  Computer  Code  Considerations 


3.3.1  Computer  code  stability 

Several  factors  must  be  considered  when  designing  and  using  a 
finite-difference  computer  code.  The  time-step  size  is  the  most 
crucial  of  these  factors.  The  time  step  must  be  smaller  than  the 
smallest  time  "response"  of  any  element  in  the  model  as  defined  in 
Eq.  3.28  from  Dusinberre  [8,  p.  13].  Equations  3.29,  3.30,  and  3.31 
show  this  response  time  for  an  external  (closest  to  the  liquid  metal 
flow)  PCM  element,  an  element  at  the  center  of  the  PCM  tube,  and  a 
liquid  metal  region  element,  respectively. 
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Due  to  the  higher  conductivity  of  the  liquid  metal,  the  time  step 
size  for  this  region  was  always  several  orders  of  magnitude  smaller 
than  that  needed  for  the  PCM  region.  This  portion  of  the  code  was 
therefore  rewritten  in  an  implicit  manner.  The  stability  of  the 
implicit  formulation  is  not  as  easily  defined  as  that  of  the  explicit 
formulation.  The  implicit  approach  is  expected  to  be  stable  at  very 
large  time  steps  but  the  only  way  to  test  the  stability  is  to  vary  the 
time-step  size  and  examine  the  effect  on  the  results.  Therefore,  the 
code  was  run  with  six  radial  subdividions  and  a  variety  of  time-step 
sizes.  The  maximum  time-step  size  for  this  examination  was  set  by  the 
requirements  of  the  explicit  portion  of  the  code.  The  results  of  this 
test  are  shown  in  Fig.  3.3  in  the  form  of  the  liquid  metal  temperature 
profile  after  100  s.  These  consistent  shape  of  these  curves  indicates 
that  the  solution  is  indeed  stable  for  the  time  step  sizes  under  con¬ 
sideration.  The  accuracy  of  the  liquid  metal  exit  temperature  is 
acceptable  considering  other  approximations  used  in  the  finite- 
difference  model. 

Radial  and  axial  step  sizes  were  examined  in  a  similar  manner. 

The  number  of  radial  subdivisions  was  varied  from  3  to  24.  The  axial 
step  size  was  automatically  adjusted  so  that  the  perimeter  and  cross- 
sectional  area  for  a  PCM  element  were  of  the  same  order  of  magni¬ 
tude.  The  finite  difference  code  automatically  stops  when  the  liquid 
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Fig.  3.3.  Influence  of  Time-Step  Size  on  the  Stability  of  the 
Implicit  Solution. 


metal  exit  temperature  exceeds  1100°K.  The  time  at  which  this  condi¬ 
tion  was  met  was  very  sensitive  to  the  number  of  radial  subdivisions, 
varying  from  964  seconds  at  3  subdivisions  down  to  697  seconds  at  24 
subdivisions.  The  results  shown  in  Figs.  3.4  and  3.5  reflect  the  PCM 
condition  at  the  stopping  point  determined  by  the  liquid  metal  exit 
temperature  limit. 

Figure  3.4  shows  the  effect  of  element  size  on  the  dimensionless 
PCM  temperature  at  the  outer  surface  of  the  tube  and  in  the  center  of 
the  tube.  The  outer  surface  temperature  is  relatively  unaffected  by 
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Fig.  3.4.  Effect  of  Element  Size  on  PCM  Dimensionless  Tempera¬ 
ture  Profile. 

the  number  of  divisions  but  the  temperature  at  the  center  of  the  tube 
is  strongly  affected.  Figure  3.5  shows  the  effect  of  element  size 
on  0*,  the  dimensionless  energy  storage  of  the  tube  and  the  Fourier 
number  corresponding  to  the  time  at  which  the  liquid  metal  exit  tem¬ 
perature  exceeds  1100°K.  Based  on  an  examination  of  these  figures, 
nine  subdivisions  were  used  for  all  subsequent  tests.  This  is  a  com¬ 
promise  between  accuracy  and  efficiency.  The  computer  code  running 
time  increases  at  a  rate  roughly  proportional  to  the  number  of  radial 
divisions,  from  14  Cray  CPU  minutes  for  M  ■  3  to  126  Cray  CPU  minutes 
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Fig.  3.5.  Effect  of  Element  Size  on  Test  Duration  (F  )  and 
Energy  Storage  (0*). 

for  M  a  24.  At  nine  subdivisions,  the  predicted  energy  storage  is 
within  13%  and  time  value  is  within  10%  of  the  more  accurate  tests. 
This  accuracy  is  therefore  considered  acceptable  considering  the  wide 
range  of  potential  design  parameters  discussed  in  Sect.  1.2.  The 
final  computer  runs  for  the  enhanced  configurations  used  from  53  to 
153  CPU  minutes  on  the  Cray  computer. 


> 
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I 

|  3.3.2  Computer  code  development 

The  overall  computer  code  flow  chart  is  shown  in  Fig.  3.6.  This 
figure  and  the  results  presented  in  this  report  reflect  the  final  form 
j  of  the  computer  code.  However,  several  subjects  explored  during  the 

development  of  this  code  are  worth  discussing.  The  change  from  an  ex- 
plicit  to  an  implicit  treatment  of  the  liquid  metal  region  has  already 
been  discussed.  Another  important  facet  of  the  development  included 
the  addition,  and  ultimately  the  deletion,  of  a  segment  of  the  code 
that  modeled  the  PCM  containment  tube  wall. 

The  inclusion  of  the  tube  wall  complicated  the  code  in  several 
ways.  The  code  was  found  to  be  better-behaved  if  the  cross-sectional 
areas  of  each  finite  difference  node  were  of  the  same  order  of  magni¬ 
tude  as  the  axial  areas.  Because  the  tube  wall  was  so  much  thinner 
than  the  PCM  region,  achieving  such  a  balance  was  very  difficult. 
Additionally,  at  the  juncture  of  the  PCM  and  the  wall,  it  was  neces¬ 
sary  to  know  the  PCM's  conductivity,  which  in  turn  required  knowledge 
of  the  phase  at  that  point  in  space  and  time.  This  somewhat  nullified 
the  advantages  inherent  in  the  Kirchoff  temperature  approach.  Addi¬ 
tionally,  the  tube  wall's  high  conductivity  and  small  physical  dimen¬ 
sions  led  to  very  small  time  steps,  on  the  order  of  those  required  for 
the  liquid  metal  region.  Moreover,  after  going  through  all  these 
correction  procedures  to  the  extent  possible,  the  temperature  differ¬ 
ence  across  the  tube  wall  was  always  much  smaller  than  the  difference 


between  the  liquid  metal  and  PCM  temperatures 
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Fig.  3.6.  Flow  Chart  for  Finite  Difference  Computer  Code. 


After  deciding  to  remove  the  tube  wall  portion  of  the  model  from 
the  code,  several  caveats  were  added  to  the  code.  The  code  now  calcu¬ 
lates  the  thermal  storage  capacity  of  the  tube  wall  and  compares  it  to 
the  thermal  storage  capacity  of  the  PCM.  This  value  is  printed  out 
and  an  excessively  large  value  will  stop  the  code's  execution.  Also, 


the  convective  heat  transfer  coefficient  is  compared  to  the  overall 
heat  transfer  coefficient  (Eq.  3.32)  including  the  wall.  This  value 
is  also  printed  out  and  can  stop  the  code's  execution.  These  precau¬ 
tions  will  be  of  value  if  other  tube  wall  thicknesses  or  materials  are 
examined  at  a  future  time. 


h 


overall 


(Eq.  3.32) 


3.4  Enhanced  Configurations 


As  will  be  demonstrated  in  Sects.  4.1  and  4.2,  the  limiting 
factor  of  the  PCM  energy  storage  system  is  the  poor  conductivity  of 
the  PCM  material.  This  conductivity,  as  was  shown  in  Table  1.1,  is 
about  one-tenth  the  conductivity  of  either  the  liquid  metal  or  the 
tube  wall.  If  the  PCM  conductivity  were  higher,  the  tube  radius  could 
be  larger,  the  number  of  tubes  could  be  smaller,  and  the  liquid  metal 
flow  rate  could  likely  be  lower. 

Three  enhancement  methods  are  considered  in  this  study.  They 
involve  replacing  some  portion,  measured  as  a  volume  percent,  of  the 
PCM  with:  1)  a  reticulated  metal  insert,  2)  thin  fins,  or  3)  liquid 
lithium  (to  make  a  metal  salt/molten  metal  slurry).  A  mixed  enthalpy 
approach  is  used  to  model  these  enhancements  with  the  finite  differ¬ 
ence  code. 


The  mixed  enthalpy  approach  is  based  on  an  assumption  that  the 
enhanced  PCM  mixture  behaves  as  a  homogenous  material.  Knowles  and 
Webb  (9]  have  defined  optimum  fin  dispersion  requirements  for  which 
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this  assumption  is  rigorously  true.  These  dispersion  requirements  are 
based  on  two  conditions:  1)  the  fins  must  be  perpendicular  to  the 
heat  transfer  face  and  2)  the  fins  must  be  thin  enough  that  the 
thermal  resistance  along  the  fins  must  be  greater  than  the  resistance 
across  the  PCM  layers.  For  cylindrical  geometries,  this  means  that 
the  enhancement  material  should  be  arranged  perpendicular  to  the  outer 
face  of  the  cylinder  and  the  concentration  should  be  inversely  propor¬ 
tional  to  the  radius.  The  applicability  of  this  approach  for  thin 
fins  has  been  experimentally  confirmed  by  Knowles  and  Webb  [9], 

The  mixed  enthalpy  method  is  implemented  by  calculating  new 
"mixed"  values  for  the  density,  heat  capacity,  conductivity,  and 
latent  heat  of  the  PCM  mixture  (see  Eqs.  3.33-3.36).  In  Eq .  3.34,  the 
specific  heat  of  each  material  is  first  converted  to  a  volumetric  spe¬ 
cific  heat  (used  and  empirically  documented  by  Knowles  and  Webb  [9]), 
weighted  according  the  volumetric  fraction  of  the  addition  material, 
and  finally  reconverted  to  a  mass  basis  using  the  density  of  the  mix¬ 
ture.  The  latent  heat  of  the  mixture  (Eq.  3.36)  must  reflect  not  only 
the  fact  that  a  portion  of  the  material  does  not  undergo  a  phase 
change,  but  also  the  new,  greater  density  of  the  mixture. 


P  .  =  x  +  (1  -  x)  p 

mix  add  pern 


xc  p  +  (1  —  x)  c  p 
p  add  p  pem 
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-  -^2“  (  i  _  x)  L 
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mix  pem 


(Eq.  3.33) 


(Eq.  3.34) 


(Eq.  3.35) 


(Eq.  3.36) 


Both  aluminum  and  nickel  are  available  as  a  reticulated  material 


(similar  to  steel  wool)  with  void  fractions  from  90  to  95%.  Aluminum 
has  a  much  higher  conductivity  but  unfortunately  a  melting  point  below 
the  application  under  consideration.  Nickel  was  therefore  chosen  as 
the  material  for  the  reticulated  metal  insert  and  its  properties  are 
shown  in  Table  3.3.  Both  90  and  95%  void  fractions  were  modeled. 

Table  3.3.  Enhancement 
Material  Properties 

Nickel12 

k  55.0  W/m-°K 

Cp  0.444  J/g-°K 

P  8.906  x  106  g/m3 

b 

Liquid  Lithium 

k  12.22  W/m-'K 

cp  4.6  J/g-°K 

p  4.65  x  jo3  g/m3 

a  ~ 
c  and  p  at 

20°C,  kPat  600°K  [3], 

&at  973°K  [14], 

The  fins  were  made  out  of  the  same  material  as  the  tube  wall, 
304SS.  The  properties  for  this  material  were  listed  in  Table  1.1. 

Fin  volume  fractions  from  5  to  50%  were  modeled. 

The  metal/metal-salt  mixture  is  generated  by  adding  liquid 
lithium  to  the  lithium  hydride  salt.  This  addition  not  only  enhances 
the  heat  transfer  but  is  also  valuable  because  it  significantly  re¬ 
duces  the  PCM's  hydrogen  pressure  as  can  be  seen  in  Fig.  3.7.  This 
reduction,  in  turn,  simplifies  the  PCM  containment  problem.  Volume 
fractions  of  liquid  lithium  from  10  to  50%  were  modeled. 


4.  RESULTS 


4. 1  Unenhanced  System  Configurations 

Initial  system  configuration  estimates  were  based  on  the  results 
of  the  approximate  analysis  described  in  Chap.  2.  An  initial  screen¬ 
ing  of  the  candidate  systems  reported  in  Table  2.1  was  made  with  only 
three  radial  subdivisions  to  save  computer  time.  Only  one  of  these 
configurations  was  able  to  satisfy  the  thermal  load  for  the  target 
test  duration  of  1000  seconds.  This  configuration  was  the  one  with 
the  smallest  value  of  RQ  and  an  Rf/R0  ratio  of  3.  Therefore,  another 
set  of  trial  cases  was  generated  with  still  smaller  tube  radii,  recog¬ 
nizing  that  as  the  number  of  radial  subdivisions  was  raised  to  nine 
(see  Sect.  3.3.1  for  a  discussion  of  accuracy  requirements),  the  ther¬ 
mal  load  would  be  satisfied  for  a  still  shorter  period  of  time.  This 
second  set  of  cases  is  shown  in  Table  4.1.  The  number  of  tubes  re¬ 
quired  for  these  estimates  was  based  on  the  thermal  storage  capacity 
of  both  the  PCM  and  the  liquid  metal  within  the  annular  region. 

All  of  these  trial  configurations  were  able  to  satisfy  the  ther¬ 
mal  storage  demand  for  the  target  operating  period  of  900  to 
1000  seconds.  However,  the  system  energy  density  was  low  and  only  a 
small  portion  of  the  PCM  was  melted.  Therefore,  for  two  configura¬ 
tions,  the  number  of  tubes  was  then  reduced  below  the  estimated 
requirements  to  increase  the  energy  storage  load  on  each  tube.  This 
increased  both  the  energy  storage  density  and  the  proportion  of  energy 
stored  in  the  form  of  latent  heat,  as  is  shown  in  Fig.  4.1.  For  two 
of  these  tube  configurations,  the  liquid  metal  exit  temperature 
exceeded  the  allowable  exit  temperature  at  operating  periods  of  900  to 
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Fig.  4.1.  Energy  Stored  per  Tube  vs  Number  of  PCM  tubes. 


950  seconds,  as  noted  on  the  table.  The  maximum  operational  period  is 


therefore,  an  important  factor  in  defining  the  maximum  energy  storage 


density  attainable.  The  last  configuration  listed  on  Table  4.1  was 


chosen  as  a  basis  of  comparison  for  the  enhanced  configurations  in 


Sect.  4.2.  It  is  important  to  note,  however,  that  this  configuration 


does  not  represent  an  optimum,  as  is  discussed  below. 


The  length  to  diameter  ratio  had  no  effect  on  the  energy  storage 


density  for  these  configurations  because  of  the  way  the  required  num¬ 


ber  of  tubes  was  estimated.  An  increase  in  this  ratio  increases  the 
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proportion  of  energy  stored  as  latent  heat  for  those  cases  with  a 
lower  energy  storage  density.  However,  for  the  cases  where  the  ther¬ 
mal  load  per  tube  is  increased,  the  length  to  diameter  ratio  has 
almost  no  effect  on  the  proportion  of  latent  energy  stored.  This 
parameter  is  apparently  important  only  for  those  cases  where  the  ther¬ 
mal  load  is  so  small  that  the  melting  front  remains  very  close  to  the 
surface  of  the  PCM  tube  throughout  the  test  period. 

The  high  (about  60%)  proportion  of  sensible  heat  storage  in  the 
PCM  is  largely  attributable  to  the  difference  between  the  initial  sys¬ 
tem  temperature  of  700°K  and  the  PCM  melt  temperature  of  962°K. 

The  amount  of  energy  stored  in  the  increased  internal  energy  of 
the  liquid  metal  heat  transfer  fluid  was  very  high,  about  1.5  times 
the  amount  of  energy  stored  in  the  PCM  tube.  This  can  be  seen  in 
Table  4.1  where  the  system  energy  storage  and  system  energy  density 
include  the  liquid  metal  inventory  within  the  annular  region.  Con¬ 
sidering  this  ineffective  distribution  of  energy  storage,  several 
other  configurations  were  tested,  with  smaller  radii  and  with  smaller 
Rf/R0  ratios.  Most  of  these  cases  were  unable  to  satisfy  the  operat¬ 
ing  period  requirements.  However,  one  test  case  did  satisfy  this  re¬ 
quirement  and  shows  that  the  system  configurations  discussed  above  are 
far  from  optimum.  This  case  was  identical  to  the  fourth  case  listed 
in  Table  4.1  except  that  the  Rf/RQ  ratio  was  decreased  to  2.0.  The 
ratio  of  latent  to  sensible  energy  storage  in  the  PCM  tube  increased 
from  0.029  to  0.504,  9*  increased  from  0.698  to  1.246,  the  system 
energy  density  increased  from  421  to  973  J/g,  and  the  ratio  of  energy 


stored  in  the  liquid  metal  to  that  stored  in  the  PCM  dropped  from  1.5 
to  0.41. 

Based  on  this  example,  and  on  the  large  number  of  system  param¬ 
eters  that  could  be  varied,  the  system  performance  reported  in  this 
study  should  be  used  only  to  measure  the  relative  value  of  system 
enhancements,  not  as  an  absolute  measure  of  potential  system  per- 
f  ormance . 

4 . 2  Enhanced  System  Comparisons 

The  three  enhancement  options  were  described  in  Sect.  3.4  and 
include  the  insertion  of  a  reticulated  nickel  material,  the  addition 
of  thin  fins,  and  the  addition  of  liquid  lithium  to  the  PCM.  These 
enhancements  had  a  modest  effect  on  system  performance,  increasing 
both  the  length  of  time  for  which  the  thermal  load  could  be  met  and 
the  total  energy  stored.  The  decrease  in  energy  storage  density 
varies  greatly  among  the  enhancement  options.  The  following  compari¬ 
sons  are  based  on  a  maximum  operating  period  of  1000  seconds,  as 
described  in  the  problem  definition  in  Sect.  1.2.  Section  4.3  will 
describe  some  alternative  bases  for  comparison. 

The  operational  period  of  the  unenhanced  system  was  947  sec¬ 
onds.  The  addition  of  thin  fins  and  the  insertion  of  reticulated 
nickel  material  increased  the  operational  period  beyond  1000  sec¬ 
onds.  The  addition  of  liquid  lithium  increased  the  operational  period 
to  960  seconds  for  a  10%  addition;  however,  the  operational  time  for 
larger  additions  was  decreased.  Figure  4.2  summarizes  the  effect  of 
these  enhancement  options  on  the  operating  period. 
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Fig.  4.2.  Effect  of  Enhancements  on  Test  Durations. 


The  amount  of  energy  stored  was  increased  by  up  to  6%  by  the  use 
of  a  reticulated  nickel  insert  or  by  the  addition  of  thin  fins,  as 
seen  in  Fig.  4.3.  The  energy  storage  benefits  of  the  90  and  95%  void 
fraction  nickel  inserts  are  the  same.  Two  thirds  of  the  energy  stor¬ 
age  increase  attainable  by  the  addition  of  thin  fins  is  available  at 
only  a  10%  volume  addition.  The  energy  storage  doesn't  increase  at 
all  when  the  volume  fraction  of  the  fin  material  is  increased  further 
from  30  to  50%.  The  addition  of  liquid  lithium  increased  the  energy 
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Fig.  4.3.  Effect  of  Enhancements  on  Total  Energy  Storage. 


storage  a  small  amount  for  additions  of  10%  or  less.  Larger  amounts 
of  liquid  lithium  actually  decreased  the  energy  storage  capacity. 

This  decrease  is  attributable  to  the  low  density  of  liquid  lithium, 
lower  even  than  the  PCM.  Since  the  volume  of  each  tube  is  held  con¬ 
stant,  increasing  the  proportion  of  liquid  lithium  decreases  the  mass 


of  the  PCM  mixture.  The  opposite  is  true  for  the  other  enhancement 


For  the  specified  operating  period  range,  the  energy  storage 
density  is  decreased  by  the  use  of  fins  and  the  reticulated  nickel 
insert  but  was  relatively  unchanged  by  the  addition  of  liquid  lithium, 
see  Fig.  4.4.  The  thin  fin  and  nickel  additions  of  10%  decrease  the 
energy  storage  density  by  about  7%. 

Figure  4.5  shows  that  a  reticulated  nickel  insert  with  a  95%  void 
fraction  offers  the  best  compromise  between  an  increase  in  storage 
capacity  and  a  relatively  small  decrease  (~1%)  in  the  energy  storage 
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Fig.  4.4.  Effect  of  Enhancements  on  System  Energy  Density. 
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Fig.  4.5.  System  Energy  Density  vs  Total  Energy  Stored  for  En¬ 
hanced  Systems. 


density  for  the  specified  operating  period  range.  An  addition  of  10% 
liquid  lithium  increases  both  the  energy  storage  density  and  the  total 
energy  storage  capacity,  but  only  by  a  very  small  amount.  However, 
this  option  has  the  added  benefit  of  reducing  the  hydrogen  pressure  of 
the  lithium  hydride  PCM,  as  was  discussed  in  Sect.  3.4. 

The  ratio  of  energy  stored  in  the  form  of  latent  heat  to  the 
sensible  heat  within  the  PCM  varies  from  0.43  in  the  unenhanced  con¬ 
figuration  to  0.56  in  the  90%  void  fraction  reticulated  nickel  insert 


configuration,  as  Table  4.2  shows.  If  all  the  unenhanced  PCM,  from  an 
initial  temperature  of  700°K,  were  melted  and  at  the  melting  tempera¬ 
ture,  this  ratio  would  equal  a  value  of  1.3.  However,  much  of  the  PCM 
remains  unmelted  at  the  end  of  the  operating  period.  Also,  the  PCM 
nearest  the  surface  heats  up  above  the  melting  temperature,  further 
increasing  the  proportion  of  sensible  heat  storage. 

The  position  of  the  phase  change  front  at  the  end  of  the  operat¬ 
ing  period  is  shown  in  Figs.  4.6  -  4.8.  In  these  figures,  the  PCM 
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Fig.  4.6.  Phase-Change  Front  Positions  for  Reticulated  Nickel 
Enhancements . 


/•ij 


M 


•ft, 

•*$ 


lo. 

1 

v.m: 


I1!.**! 


SSI 

I 

I 

$ 

r 

I 


ijjc] 

i 

S5S 


& 

iyc| 

»s 


>% 

bO 

w*  - 

g  =  2 

•  „  B 

4)  _ 
O  U.  ^ 

o  -d 

O  u  “ 
a)  5 


*  3 
5  So" 

Q  W  C 
B  0)  *H 

„  c 

2  O)  T3 

3  « 

i-1  o 


<T  <T  *3 


O  ^ 

•m  S  «  **  y 
°  ^  r, 
o -5 ? 8,* 
r  S  g  g  2  « 
<3  «  S  "  2  - 

tn 


0) 

r_,<V*a  ^ 
«  68  £." 
o  a  P 
H  s »- 


O  O 
sr  in  in 


N  O  PI  CO 
^  m  in  4 


o  o  o  o 


"3  A) 
’XJ  bO  45io  s 

r  u  ^  o  -e 

u  a>  o  -«  2 

ft-  C  4J  V-*' 

<u  03  zr 


3  -v.  ^ 

e  bo  W  ^ 
a>  Jt-h  bo 

«2c*-) 

C/5  W  *3 


^  O' 
g  00  <r 

«  v©  v© 


Os  v©  — » 
rx.  r^>  r^ 


n  oo  tn 
CM  m  O 


— 4  CO  CO 

r*^  r*«. 


CM  CM  <* 

O  O'  P' 
r^.  v©  sO 


o  n  <r 
v©  <r  — « 

3v  Cs  CT' 


m  in  n  n 
00  00  00  00 


>©  O0  Id  N 

n  tn  oo  m 


^  oo  >©  N 


N  in  (Tv  -H 

r*»  r**  r*.  oo 


s©  cm  oo  ^ 
m  mt 

v©  -©  in  ^ 


00  CM  O  O 
r*»  O'  O  O 
O'  O'  o  o 


•  45 

O  X 

O  U 

O  <0 

sf  45 

u 

N 

<u 

CO  I- 

CO  3 

ft-  AJ 

M  <0 

ft-  u  • 

55  45  »H 

a 

•  0  <0 

O  45  J* 

m  4ui 

■  •*-*  ^ 

•H  3 

Q  X  *J 

^  4) 

►J  t3 

»-»  C 

•  W  51 

00  4J 


5  C  H 

S  O  CO 


CO  >  CM 

e  s 

u  cn 
a>  3io 
X  y-*  l 
H  •*-»  O 


Mf 

*3 

bO 

V- 

O 

■H 

4/ 

45 

• 

3 

U 

4J 

O 

cr 

3 

•H 

H 

r^. 

cjn 

V© 

crv 

O 

V© 

N© 

r-. 

1 

H 

05 

CM 

<f  v© 

m 

in 

o 

cn 

ON 

ON 

«J 

£ 

• 

•  • 

• 

• 

• 

• 

• 

• 

• 

> 

c 

45 

U 

o 

o  o 

o 

o 

o 

o 

o 

4 ) 

e 

ft- 

•> 

x: 

cn 

5 

T3 

•3 

•H 

3 

l 

05 

3 

O 

CL 

O' 

C 

o 

o 

•H 

fl 

02 

U 

H 

x: 

05 

c 

O' 

m  o 

CM 

cn 

<r 

05 

45 

• 

«  • 

• 

• 

• 

• 

• 

• 

■ 

02 

«-M 

45 

o 

— •  CM 

*-4 

w- 4 

CM 

•H 

CM 

<n 

05 

•3 

» 

TJ 

3 

B 

O 

»M 

c 

fi 

3 

o 

in 

c 

O 

u 

•H 

0) 

• 

4) 

3 

o 

B 

•3 

3 

45 

3 

1 

U 

0) 

H 

c; 

o 

O 

c 

02 

u 

CO 

»-« 

o 

in  O 

o 

o 

O 

m 

o 

o 

O 

<3  - 

-Q 

3 

*X3 

cn 

m 

cn 

m 

mm 


Includes  energy  stored  In  the  form  of  Increased  Internal  energy  of  liquid  metal  in  the  annular  region 
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Fig.  4.7.  Phase-Change  Front  Positions  for  Liquid  Lithium  En¬ 
hancements. 

above  the  line  is  melted  and  the  PCM  below  the  line  is  solid.  The 
stair-step  appearance  of  these  plots  is  caused  by  the  nature  of  the 
finite  difference  solution  and  the  size  of  the  radial  divisions.  The 
front  for  most  of  the  cases  only  penetrates  about  25%  of  the  way  into 
the  PCM.  However,  at  an  r/RQ  value  of  .75,  about  44%  of  the  PCM 
volume  is  in  the  melted  zone.  The  center  of  the  PCM  is  melted  only 
for  very  large  additions  of  enhancement  materials.  Such  large  addi¬ 
tions  are  not  supported  by  the  energy  storage  and  energy  density  con¬ 
siderations  already  discussed. 
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Fig.  4.8.  Phase-Change  Front  Positions  for  304SS  Thin  Fin  En¬ 
hancements. 


It  is  interesting  to  note,  however,  that  doubling  the  volume  per¬ 
cent  of  enhancement  material  from  5  to  10%  has  a  relatively  small 
effect  on  the  interface  location  for  both  the  reticulated  nickel  and 
thin  fin  cases,  as  does  tripling  the  enhancement  material  from  10  to 
30%  in  the  liquid  lithium  case.  However,  tripling  the  fin  volume  from 
10  to  30%  has  a  very  large  effect,  moving  the  interface  all  the  way  to 
the  center  of  the  PCM  tube  for  z/L  <0.2.  Increasing  the  fin  volume  by 
another  67%  causes  the  melting  front  to  reach  the  center  for  80%  of 


the  length  of  the  tube.  Also,  Increasing  the  liquid  lithium  fraction 
from  30  to  50%  (only  a  67%  increase)  causes  a  tremendous  increase  in 
melt  front  propagation  for  z/L  <0.1. 


Figure  4.9  shows  the  dimensionless  temperature  response  of  the 
central  (r=0,  z/L=.5)  PCM  node  vs.  the  dimensionless  time  variable 
which  reflects  the  varying  Stefan  numbers  of  the  enhanced  PCM. 

The  nondimens ional  temperature  response  is  slightly  slowed  by  addi¬ 
tional  amounts  of  enhancement  material  (i.e.,  3<Ji/3FoSte  is  decreased 
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Fig.  4.9.  Nondimensional  Time  Response  of  Central  PCM  Node  for 
Enhanced  Configurations. 


as  the  thermal  diffusivity  is  increased).  The  nondimensional  behavior 


reflects  the  fact  that  the  FoSte  time  parameter  is  proportional  to  the 
increasing  thermal  conductivity  of  the  enhanced  materials.  This 
increasing  thermal  conductivity  stretches  the  time  scale,  causing  this 
counter-intuitive  trend.  However,  for  the  real  time  data  shown  in  Figs. 
4.10-4.12,  the  time  response  is  improved  (i.e.  3T/3t  is  increased)  by 
the  addition  of  the  enhancement  materials. 
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Fig.  4.10.  Temperature  Response  of  Central  PCM  Node  for  Systems 
Enhanced  With  Reticulated  Nickel. 


4.3  Alternative  Bases  for  Comparison 


The  performance  comparisons  reported  in  Sect.  4.2  were  based  on  a 
strict  adherance  to  the  problem  definition.  Specifically,  operating 
periods  longer  than  1000  seconds  were  not  considered  even  though  some 
of  the  enhanced  systems  were  capable  of  receiving  and  storing  further 
thermal  energy  without  exceeding  the  liquid  metal  inlet  and  outlet 
temperature  constraints.  If  the  application  of  this  thermal  storage 
system  was  well  defined  and  unlikely  to  change,  this  basis  for  compar¬ 
ison  would  be  unquestioned.  However,  given  the  vagueness  of  the  cur¬ 
rent  applications  requirements,  further  discussion  of  the  performance 
of  these  systems  is  in  order. 

If  the  comparison  were  based  solely  on  the  system's  physical 
dimensions  and  on  the  inlet  and  outlet  liquid  metal  temperature  con¬ 
straints,  several  conclusions  reached  in  Sect.  4.2  would  be  different 
as  can  be  seen  in  Table  4.3. 

For  the  reticulated  nickel  insert,  the  most  important  change  is 
that  now  the  total  energy  stored  for  the  reticulated  nickel  insert 
continues  to  increase  from  187  MJ  to  197  MJ  as  the  volume  percent  of 
nickel  is  increased  from  5  to  10%.  As  the  void  fraction  drops  from  95 
to  90%,  the  system  energy  density  is  decreased  by  only  20  J/g  rather 
than  the  sharper  drop  of  almost  40  J/g  with  the  1000  second  limit. 

This  would  change  the  conclusion  that  the  95%  void  fraction  material 
is  clearly  more  effective  than  the  90%  void  fraction  material.  The 
choice  between  these  two  void  fractions  then  becomes  a  more  subtle 
trade-off  between  increased  storage  capacity,  lengthened  operating 
period,  and  slightly  decreased  system  energy  density. 


The  effect  of  this  approach  on  the  fin  results  is  not  as  dramatic 
because  the  overall  system  energy  density  is  still  decreased  substan¬ 
tially  by  the  use  of  such  fins.  However,  it  becomes  apparent  that  a 
maximum  energy  storage  capacity  occurs  somewhere  between  10  and  50% 
fin  volume.  The  operating  period  is  also  slightly  longer  within  this 
range . 

Considering  the  changeable  nature  of  these  conclusions,  some  word 
of  caution  should  accompany  the  results  of  this  study.  It  would 
appear  that  the  most  appropriate  system  comparison  should  be  based  on 
optimized  parameters  for  each  enhanced  system.  However,  such  an 
optimization  strategy  is  beyond  the  scope  of  this  study.  The  develop¬ 
ment  of  such  a  strategy  should  be  pursued  if  this  type  of  thermal 
storage  system  remains  a  viable  candidate  in  the  overall  system 


5.  CONCLUSIONS  AND  RECOMMENDATIONS 


An  approximate  model  of  the  PCM  energy  storage  system  was  used  to 
identify  preliminary  physical  dimensions  for  the  specified  applica¬ 
tion.  This  specified  application  included  a  short  cnarging  period  and 
the  use  of  a  liquid  metal  heat  transfer  fluid.  A  finite  difference 
computer  code  was  developed  and  used  to  evaluate  the  effect  of  three 
heat  transfer  enhancements  on  the  PCM  energy  storage  system. 

For  a  maximum  operating  period  of  1000  S,  the  use  of  a  95%  void 
fraction  reticulated  nickel  insert  was  found  to  increase  the  storage 
capacity  of  the  system  by  ~6%  with  a  decrease  in  the  system  energy 
density  of  only  ~1%.  For  the  given  physical  system,  the  storage 
capacity  can  be  increased  still  further  if  the  operating  period  is  not 
constrained.  The  addition  of  10%  liquid  lithium  was  found  to  cause 
minor  increases  in  both  storage  density  and  storage  capacity  with  the 
added  benefit  of  reducing  the  hydrogen  pressure  of  the  lithium  hydride 
PCM. 

The  unenhanced  system  that  formed  the  basis  of  comparison  for  the 
various  enhancements  was  not  an  optimized  configuration.  Based  on  a 
limited  examination  of  other  potential  system  configurations,  signif¬ 
icant  performance  improvements  are  possible.  Such  an  optimization 
effort  should  be  pursued  if  these  configurations  continue  to  be  con¬ 
sidered  as  candidate  solutions.  It  also  appears  that  a  combination  of 
a  small  addition  of  the  liquid  lithium  might  be  considered  in  conjunc¬ 
tion  with  the  reticulated  nickel,  thereby  gaining  both  the  increased 
storage  capacity  and  the  reduced  hydrogen  pressure  benefits  of  these 


two  enhancements. 
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